Fragility and elasticity - a new perspective 
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The fragility (the abnormally strong temperature dependence of the viscosity) of highly viscous 
liquids is shown to have two sources. The first is the temperature dependence of the barriers between 
inherent states considered earlier. The second is the recently discovered asymmetry between the 
actual inherent state and its neighbors. One needs both terms for a quantitative description. 

PACS numbers: 64.70.Pf, 77.22. Gm 



Though there is as yet no generally accepted expla- 
nation of the flow in highly viscous liquidaii^i^, it seems 
clear that its description requires the passage of high en- 
ergy barriers between inherent states, i.e. local structural 
minima of the potential energy- . According to the elastic 
models^, the fragility stems from a proportionality of the 
height V of these barriers to the infinite frequency shear 
modulus G, which in the highly viscous liquid decreases 
strongly with temperature. 

The present paper shows that this is only part of the 
truth. There is a second source of fragility in the newly 
discovered^ asymmetry of about AksT between the ac- 
tual inherent state and its neighbors, possibly due^ to the 
elastic distortion accompanying a structural rearrange- 
ment (the "Eshelby backstress"— ) . As will be seen, the 
quantitative explanation of the fragility of six different 
glass formers requires just this specific explanation of the 
asymmetry. 

The usual measure of the fragility of a glass former is 
the logarithmic slope of the relaxation time Tq of the flow 
process 



m^d\ogT^/d{Tg/T)\T^ 



(1) 



where the glass temperature Tg is defined as the temper- 
ature with Tq = 1000 s. 

It is useful to relate Tq to a critical barrier Vc via the 
Arrhenius relation 

r„ =roexp(14ABT), (2) 

where the microscopic attempt frequency is at 10""'^^ s, 
sixteen decades faster than the flow process at the glass 
temperature. The fragility index I is defined^, by the 
logarithmic derivative / = — din Vc/rflnT, taken at Tg. 
Then 



16(7 + 1) 



(3) 



where the factor refiects the sixteen decades between mi- 
croscopic and macroscopic time scales. / is a better mea- 
sure of the fragility than m, because it does not contain 
the trivial temperature dependence of any thermally ac- 
tivated process. 

The elastic models^ postulate a proportionality be- 
tween the flow barrier Vc and the infinite frequency shear 



modulus G. One can again define a dimensionless mea- 
sure r for the temperature dependence of G in terms 
of the logarithmic derivative F = — dlnG/dlnT at Tg. 
Then the elastic models2, postulate / = F. 

In order to check this relation, one needs measure- 
ments of both quantities. The flow relaxation time 
is relatively easy to measure, but the determination of 
the high frequency shear modulus is by no means triv- 
ial. It requires the measurement of the density p and 
the high frequency transverse sound velocity vt. Con- 
sequently, the logarithmic derivative F is a sum of two 
terms, a larger one from the sound velocity and a smaller 
one from the density. 

In a liquid, well-defined transverse sound waves do only 
exist at frequencies which are markedly higher than the 
inverse I/t^ of the flow relaxation time. With increasing 
temperature, gets very rapidly shorter. Therefore the 
measurement of the transverse sound velocity by Bril- 
louin scattering is limited to a small temperature region 
above Tg. This, together with the poor visibility of the 
transverse Brillouin line, leads to a large error bar in the 
determination of F, usually about 20 %. 
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0.3 


Vit-4 


627 


34 
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0.88 


0.56 










glycerol 


187 


4.5 


53 


2.31 


1.0 


1.59 


0.12 


2.2 


1.2 


PB20 


173 


1.8 


84 


4.25 


2.0 


1.05 


0.12 


4.8 


2.9 


CKN 


343 


4.9 


93 


4.81 


2.6 










BPA-PC 


418 


0.8 


132 


7.25 


4.4 


1.51 


0.23 


7.8 


3.4 


PS 


375 


1.0 


138 


7.63 


4.0 


0.82 


0.21 


9.8 


5.8 


PMMA 


379 


1.9 


145 


8.06 


2.1 


0.30 


0.23 


11.7 


9.6 



TABLE I: Measured and calculated fragilities for eight glass 
formers. Vit-4 is a bulk metallic glass, PB20 is a 20:80 
mixture of 1,2-polybutadiene and 1,4-polybutadiene, CKN is 
K3Ca2(N03)7, BPA-PC is bisphenol-A-polycarbonate, PS is 
polystyrene and PMMA is polymethylmethacrylate. Refer- 
ences see refi^. 



In spite of these difficulties, six apparently reliable 
measurements of F by Brillouin scattering exist in the 
literature. The six substances are listed in Table I, to- 
gether with the bulk metallic glass Vit-4 and polystyrene, 
where F was determined from ultrasonic measurements. 
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Table I compares fragility indices / calculated from 
the temperatm'e dependence of Tq, with F. One finds 
that / is always larger than F, in several cases clearly 
beyond the estimated error bar of 20 %. This is very 
surprising, because one would have expected the op- 
posite result. The barriers of the energy landscape 
in molecular glass formers are frequently intramolecu- 
lar barriersii, much less temperature dependent than 
the van-der-Waals dominated shear modulus. The same 
holds for polymers, where the torsional barrier a^^'^'^ are 
practically temperature-independent. Thus one would 
expect / < F, but one finds F < /. Though there is a 
clear tendency of a fragility increase with increasing F, 
the postulate 14 G is obviously not sufficient to explain 
the full observed fragility. 

A second possible source of fragility is the recently dis- 
covered asymmetry between the actual inherent state and 
its neighbors^. The strength of an asymmetric relaxation 
increases with increasing temperature, because the ther- 
mal population of the upper level increases. 

In order to quantify this influence within the asymme- 
try model^, an extension of the coupling modeP, consider 
its basic definition of the characteristic multi-minimum 
parameter 



In 



N3 



Gv 



(4) 



In order to determine fc, one needs an asymmetry 
model fit of dynamical shear data for the given substance, 
preferably close to the glass temperature, but still in the 
equilibrium liquid. For this purpose, one can follow the 
recipe given in the asymmetry paper—, describing the dy- 
namical shear data with the three parameters G, Vc and 7 
(7 describes the exponential rise fo{V) ~ exp{'yV/kBT) 
of fo{V) at Vc). For G, one takes the measured high- 
frequency value from Table I. If necessary, one can add a 
gaussian to describe an eventual Johari-Goldstein peaki^, 
which requires three more parameters, height, position 
and width. In this way, fcVc and consequently I2 were 
obtained for the six glass formers in Table I for which 
dynamical shear data close to Tg exist in the literature^. 
Fig. 1 (a) and (b) show the fits of the dynamical shear 
data. 

Table I compares I2 and the supposedly full fragility 
Ii = T + I2 to the measured value /. It is immediately 
clear that one needs I2 to understand the full measured 
fragility. But sometimes Ii is a bit too large, particularly 
in polystyrene and in PMMA. This is not unexpected, be- 
cause the elastic modeP expectation V ^ G holds only in 
substances where the contribution of intramolecular bar- 
riers is negligibl&ii. In polymers, the torsional barriers do 
play a rol e^^'^'^ and do not share the temperature depen- 
dence of G. In PMMA, where the relaxation is dominated 



where T is the temperature and v is the atomic or molec- 
ular volume, cjv is a temperature-independent measure 
of the density of stable states for N atoms or molecules 
in distortion space, assumed to be constant. N must 
be large enough to meet the condition /jv = 1 for the 
breakdown of the shear modulus at Tg. The barrier Vc is 
the lowest barrier for all possible A'' to reach /at = 1. In 
terms of the definitions of the coupHng model^, the jumps 
into neighboring inherent states with barriers below Vc 
are the primitive relaxations. The asymmetry model^ 
postulates that their elastic interaction brings the shear 
modulus down to zero. 

If this is indeed so, the contribution to the fragility de- 
pends crucially on the barrier density of primitive relax- 
ations at Vc. In the model, the primitive relaxation den- 
sity is characterized by a barrier density function foiV) 
and Vc is given by the 1/3-rule 



foiV) = \. 
-J 



(5) 



According to its definition in terms of /jv, foiV) 
increases with increasing temperature proportional to 
{T/GvY. This shifts Vc downwards and so provides an 
additional fragility index I2 



1 + F - avTg 

IcVc 



(6) 



Here ay is the volume expansion coefficient above the 
glass temperature Tg and fc = /o(K)- 
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■ glycerol 192.5 K 
□ PB20 178 K 
A PS 380 K 




1 2 3 
log(co/s^^) 



FIG. 1: Fits of dynamical shear data in terms of the asym- 
metry model^ (continuous hnes) for (a) silica, BPA-PC and 
PMMA (b) glycerol, PB20 and polystyrene. Note that a high 
/c means a high value of G" /G at the a-peak. Position and 
width of the secondary peaks of BPA-PC and PMMA were 
taken from mechanical relaxation data in the glass phase. 



3 



by a side group relaxation, one should rather compare / 
with J2. In fact, in this case I2 agrees with I within the 
error bars of the three measurements involved^. 

From eq. ([S]), one sees that the second contribution to 
the fragility is weak whenever fc is high. This is often the 
case in molecular glass formers (glycerol is an exception) . 
For instance, in triphenylethylene^ /c = 5.31, markedly 
higher than the one in silica. This explains why a reason- 
able agreement with the elastic model postulate Vc ^ G 
alone has been found in molecular glassesl^iiM. But even 
in these two papers, there are several cases which clearly 
have a higher fragility. 

The elastic model postulate is the most obvious con- 
nection between fragility and fast vibrations, but it is by 
no means the only proposition in this direction. Other 
proposals relate the fragility to the nonergodicity factor 
measured in high-resolution x-ray scatteringi^ or to the 
Poisson ratio The relation to the nonergodicity fac- 
tor is understandable, because a low nonergodicity factor 
means a low level of density fluctuations, which in turn 



means that one is close to the ideal glass of the Kauzmann 
paradoxon and expects a high thermodynamic fragilityi^. 
Very recentlj^, it has been pointed out that there are 
exceptions from the nonergodicity rule due to a strong in- 
fluence of secondary relaxations^^, a reasoning which is 
parallel to the one in the present work, a second fragility 
influence which requires not only a knowledge of the fast 
motion, but also of the relaxations themselves. 

To summarize: A full quantitative understanding of 
the fragility requires the consideration of both the tem- 
perature dependence of the barriers and the influence of 
the recently discovered^ asymmetry of an actually occu- 
pied inherent state with respect to its neighbors. The 
latter plays a minor role if the density of relaxations at 
the critical flow barrier is high, but can become dominant 
if it happens to be low. The quantitative agreement sup- 
ports the validity of the explanation of the asymmetry 
in terms of an elastic distortion^, an Eshelby backstress^ 
which tends to stabilize the occupied inherent states. 
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